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The time-evolution and space-distribution of internal electromagnetic fields in heavy-ion reactions at beam 
energies between 200 and 2000 MeV/nucleon are studied within an Isospin-dependent Boltzmann-Uhling- 
Uhlenbeck transport model IBUUl 1. While the magnetic field can reach about 7 x lO'^ G which is significantly 
higher than the estimated surface magnetic field (~ lO'^ G) of magnetars, it has almost no effect on nucleon 
observables as the Lorentz force is normally much weaker than the nuclear force. Very interestingly, however, 
the magnetic field generated by the projectile-like (target-like) spectator has a strong focusing/diverging effect 
on positive/negative pions at forward (backward) rapidities. Consequently, the differential n^ JTl^ ratio as a 
function of rapidity is significantly altered by the magnetic field while the total multiplicities of both positive 
and negative pions remain about the same. At beam energies above about 1 GeV/nucleon, while the integrated 
f^ ' ratio of total n^ to ;r+ multiplicities is not, the differential n^ JTl^ ratio is sensitive to the density dependence 

^SJ . of nuclear symmetry energy iisyn,(p). Our findings suggest that magnetic effects should be carefully considered 

in future studies of using the differential n^ JK^ ratio as a probe of the £'sym(p) at supra-saturation densities. 

PACS numbers: 41.20.-q, 25.70.-z, 21.65.Ef 



I. INTRODUCTION 

O ■ Magnetic fields exist everywhere in the Universe. To set the scale and appreciate the strong magnetic fields created during 

3 , heavy-ion collisions, we first recall the magnitudes of several typical magnetic fields from various sources. Many spiral galaxies 

I '•^1 ' have magnetic fields with a typical strength of ^ 3 x 10^^ G [1] and it is estimated that the intergalactic magnetic fields presently 

\ have an intensity of about < 10^^ G [2]. Some people believe that the present magnetic field of the Universe is amplified from a 

T— I ■ seed about 10"^" G by the dynamo mechanism [3, 4] while magnetic fields up to 10^^ G might appear in the early Universe [4]. 

K*" ' The strongest magnetic field of about lO'^ G near the surfaces of magnetars [5, 6] or even higher (lO'^-lO'^ G) associated with 

^^ . the cosmological gamma-ray bursts [7] have been found from astrophysical observations. Due to the limit of tensile strength 

of terrestrial materials, the strongest man-made steady magnetic field is only about 4.5 x 10^ G. To our best knowledge, it was 

-^ , first pointed out by Rafelski and Miiller that, in addition to strong electrical fields, unusually strong magnetic fields are also 

• • created in heavy-ions collisions (HICs). In sub-Coulomb barrier Uh-U collisions, the magnetic field was estimated to be on the 

f]^ \ order of lO'^ G [8]. More recently, it has been shown by Kharzeev et al. that HICs at RHIC and LHC can create the strongest 

magnetic field ever achieved in a terrestrial laboratory [10]. For example, in noncentral AuH-Au collisions at 100 GeV/nucleon, 

the maximal magnetic field can reach about lO'^ G [10, 14]. It thus provides a unique environment to investigate the Quantum 

Chromodynamics (QCD) at the limit of high magnetic field. Indeed, the study of quark-gluon-plasma under strong magnetic 

field has attracted much attention by the high energy heavy-ion community, see, e.g., ref. [9] and references therein. In particular, 

it has been shown theoretically that [10-13] QCD topological effects in the presence of very intense electromagnetic fields, i.e., 

\^ • the "Chiral Magnetic Effect", may be an evidence of local parity violation in strong interactions. Experimentally, interesting 

5^ ' indications have been reported, see, e.g., refs. [15,16]. 

Stimulated by the interesting findings at RHIC and realizing that all transport model studies of magnetic effects have so far 
focused on high energy HICs [14, 17], we investigate in this work first the strength, duration and distribution of internal magnetic 
fields created in HICs at beam energies between 200 and 2000 MeV/nucleon. This is the beam energy range covered by several 
accelerators in the world. We then focus on identifying possible magnetic effects on experimental observables using an isospin- 
dependent Boltzmann-Uhling-Uhlenbeck (BUU) transport model IBUUll [18, 19]. We find that while the magnetic field can 
reach about 7 x lO'^ G in these reactions, it has almost no effect on nucleon observables as the Lorentz force is negligibly 
small compared to the nuclear force. Very interestingly, however, the magnetic field generated by the projectile-like (target-like) 
spectator moving forward (backward) in the center of mass frame has a strong focusing/diverging effect on positive/negative 
pions moving forward (backward). As a result, the differential n^ jw^ ratio as a function of rapidity is significandy altered by 
the magnetic field while the total W and 7r+ multiplicities remain about the same. 



o 






t Corresponding author. Bao-An_Li@tamu-commerce.edu 
'Electronic address: only.ouli@gmail.com 



The paper is organized as follows. In the next section, we outline how the internal electromagnetic fields in HICs are calculated 
in the IBUUll transport model. The characteristics of the electromagnetic fields and their effects on several experimental 
observables in intermediate energy HICs are then discussed in Section III. Finally, a summary is given at the end. 

II. THE MODEL 

In the presence of electrical and magnetic fields E and B, the BUU equation can be written as 



d P 



f{r,p,t)^I{r,p,t) (1) 



where I{r,p,t) is the collision integral simulated by using the Monte Carlo method. The electrical field E (Coulomb field) 
has already been considered in most transport models. To include consistently both the electrical and magnetic fields satisfying 
Maxwell's equations, the Lienard-Wiechert potentials at a position r and time t are evaluated according to 
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^E^r,) = 4^i:Z«p-^-^(ci?«-/^«^«) (2) 
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where Z„ is the charge number of the nth particle. R„ = r — r'„ is the relative position of the field point r with respect to the 
position r', of particle « moving with velocity v„ at the retarded time tm = f — |r — r'j(f„i)l/c- The summation runs over all 
charged particles in the reaction system. In non-relativistic cases, i.e., all particles satisfy the condition v ^ c, the Eq. (2) and 
(3) reduce to the classical expressions 
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The first equation is essentially the Coulomb's law, and the latter is the Bio-Savart law for a system of moving charges. 

To take into account accurately the retardation effects, the phase space information of all nucleons before the moment t are 
required to calculate the electromagnetic fields at that moment. Some special care is thus necessary in initializing the reaction. 
In principle, the two colliding nuclei should be initialized to come from infinitely far away towards each other on their Coulomb 
trajectories. In practice, considering the need of keeping the initial nuclei stable and the computing time low, the initial distance 
between the surfaces of the two colliding nuclei is taken as 3 fm in our calculations. We make a pre-collision phase space 
history for all nucleons assuming that they are frozen in the projectile/target moving with a center of mass velocity Vp/,, i.e., 
fi = rf + Vpi, ■ t, where rj* is the initial coordinate of the nucleon. As we shall show, comparisons of our transport model 
calculations with analytical estimates for two moving charges (target and projectile) in both relativistic and non-relativistic cases 
indicate that our method of handing the pre-collision phase-space histories of all nucleons is reasonable. 

We refer the BUU code used in this study IBUUl 1. Compared to the IBUU04 [18] where the MDI (Momentum-Dependent- 
Interaction) is used [20], besides the electromagnetic fields with retardation effects, an isospin-dependent three -body force 
[21] (instead of the standard one used in the MDI, Gogny and Skyrme effective interactions) is used. Moreover, the high- 
momentum tail of the MDI isoscalar potential is readjusted to better fit the nucleon optical potential from nucleon-nucleus 
scattering experiments. Details of these modifications and their effects on experimental observables will be presented in a 
forthcoming publication [19]. In this work, we focus on the magnetic aspect of HICs at intermediate energies. Since one of 
our main motivations here is to see whether experimental observables known to be sensitive to the /Ssymlp) is affected by the 
magnetic effects, we notice here that in the IBUUll the EsymCp) is controlled by a parameter x introduced in the three-body 
part of the MDI interaction [20, 21]. By adjusting the parameters one can mimic diverse behaviors of the £'sym(p) predicted by 
various microscopic many -body theories [50] . As an example, shown in Fig. 1 are the fisym (p ) with jc = 1 , and — 1 , respectively. 
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FIG. 1: Density dependence of nuclear symmetry energy used in the IBUUll calculations. 



m. RESULTS AND DISCUSSIONS 

In this section, we first illustrate and discuss the beam energy and impact parameter dependence of the time-evolution and 
space-distribution of magnetic field. To help understand the magnetic effect in HICs, we shall also compare the Lorentz force 
with the Coulomb and nuclear forces. We then present and discuss magnetic effects on experimental observables. 



A. Characteristics of internal electromagnetic fields in lieavy-ion reactions 



Features of the internal electromagnetic fields are independent of the symmetry energy parameter x. In this subsection, unless 
otherwise specified a value of x = 1/3 is used. We take the z {x) axis as the beam (impact parameter) direction. Based on the 
formula of magnetic field strength in Eq. (3), the dominant component of the internal magnetic field is in the y axis perpendicular 
to the reaction plane {z — x). The component in the reaction plane is negligible because of the slow motions of nucleons in the 
X ox y directions especially in the early phase of the reaction. To test our approach used in calculating the electromagnetic 
fields, we first compare the magnetic field 5,(0) at the center of mass of the reaction system calculated using the full IBUUl 1 
dynamically with those obtained under some limiting conditions for idealized situations. Shown in Fig. 2 are the values of 
By{0) for AuH-Au reactions at a beam energy of 500 AMeV and an impact parameter of b=5 and 20 fm, respectively. As a 
reference, the approximate magnetic field of 10^^ G on the surfaces of magnetars is also indicated. The legend "classical" 
and "relativistic" indicate results obtained using Eq. (5) and Eq. (3), respectively. For a comparison, we have also performed 
calculations using both Eq. (5) and Eq. (3) assuming that the projectile and target are two point charges located at their individual 
centers of masses and are moving with their initial velocities only. Results of this calculation are denoted by the "kine.". Several 
interesting observations can be made. Firstly, it is seen that the By{Q) calculated with the classical and relativistic formulas are 
very close to each other, for both the kinematic and dynamical calculations, as one expects for reactions at relatively low beam 
energies. Secondly, the dynamical IBUUll results and the kinematic estimates are very close at the beginning and the end of 
the reaction, but they are very different during the reaction phase spanned by the small balls of the same color. The magnetic 
field has contributions from the projectile-like and target-like spectators as well as charged particles in the participant region. 
Contributions from the latter, however, are very weak because of the approximately isotropic nucleon momentum distribution 
there. Once the projectile and target begins overlapping, nucleon-nucleon collisions will start transferring the participants' 
longitudinal momenta into transverse directions. Thus, the By{Qi) from the IBUUl 1 is weaker than the kinematic estimate during 
the reaction phase. We notice that the magnetic field in the x and z directions are rather weak because they only come from 
charged participants which are moving essentially randomly in all possible directions. For the very peripheral reactions with 
^7=20 fm, the two nuclei do not overlap. As one expects, thus there is almost no difference between the kinematic and dynamical 
results. The above comparisons enhance our confidence in using the IBUUl 1 model to study the internal electromagnetic fields 
and their effects in HICs. In the following, we only present results calculated with the relativistic formula and the dynamical 
IBUUll model. 

The contours of the nucleon density p/po , the magnetic field strength eBy, and the electric field strength eE^ in the x — z 
plane at f=10, 20, 30 and 40 fm/c for the 500 AMeV AuH-Au collisions at an impact parameter of b=\Q fm are shown in Fig. 3. 
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FIG. 2: (color online) Time evolutions of the magnetic field strength eBy{0) at the center of mass of the reaction system for 500 AMeV Au+Au 
reactions at ^1=5 and 20 fm, respectively. The magnitude of e6y(0) for b=20 fm is multiplied by 10 for clarity. The approximate beginning 
and ending of the overlap phase between the projectile and target are indicated by the small balls (for b =20 fm there is no overlap). 
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FIG. 3: (color online) Distributions of the nucleon density p/po (upper panel), the magnetic field strength eBy (middle panel) and the electrical 
field strength eSv (lower panel) in the x — z plane at t=lO, 20, 30 and 40 fm/c for the 500 AMeV Au+Au collisions at an impact parameter of 
b=lO fm. 



We notice that both the eBy and eB^ are plotted here in unit of MeV^ which is equal to 1 .44 x lO'^ G. For discussing the spatial 
distribution of the electromagnetic fields, we can divide the space into three zones in terms of the x coordinate; the outside- 
zone where \x\ > 15 fm; the spectator-zone where 5 fm < |x| < 15 fm; and the overlap-zone where |.^| < 5 fm. As mentioned 
above, the electromagnetic fields come from both the spectators and participants. In the outside-zone, the spectator near the field 
point generates a stronger magnetic field in the negative y-direction while the other spectator farther away generates a weaker 
magnetic field in the positive y-direction. The superposition leads to a magnetic field points to the negative y-direction. On the 
other hand, the electric field eEx in the outside-zone includes contributions from all charges. Its sign is the same as the sign of 
the x-coordinate of the field point. In the overlap-zone, the magnetic fields generated by the two spectators will superimpose 
constructively since they are all in the positive y-direction, while the magnetic fields generated there by the moving charges in 
the participant region will largely cancel each other The strength of the magnetic field peaks when the two nuclei have reached 
the maximum compression. It then drops when the spectators depart from each other. The signs of the electric filed in the 
x-direction generated by the two spectators are always opposite, leading to the very weak electrical field in the participant region 
where the magnetic field is the strongest. 

Next, we explore the impact parameter and beam energy dependence of the magnetic field at the center of mass of the reaction 
system. Shown in the right panel of Fig. 4 is the impact parameter dependence of eBy{Q). The strength of magnetic field grows 



25 



^ 20 

3 

nj 
O 

^ 15 

o 



CQ 



5 - 



- 



I ' I ' 
Au+Aub=10fm 

200AMeV_ 

SOOAMev" 

700AMeV 

lOOOAMeV 

1500AMeV~ 

Magnetar 





I ' I ' 
500 AMeV Au+Au- 

b= Ofm 

b= 2fm 

b= 5fmH 

b= 8fm 

b=10fm 

b=12fm- 

b=15fra_ 
b=20 fm 
Magneta? 



10 



20 30 
t (fm/c) 



40 



50 



10 



20 30 
t (fm/c) 



40 



50 



FIG. 4: (color online) (Right panel) Impact parameter dependence of eSv(O) for the 500 AMeV Au+Au collisions. (Left panel) Beam energy 
dependence of eSy (0) for the Au+Au collisions with fc = 1 fm. Durations of the overlap between the proj ectile and target are indicated roughly 
by the open circles (for b =15, 20 fm there is no overlapping). The thin black dashed line is for the magnetic field strength at the surfaces of 
magnetars. 



with increasing impact parameter b up to about b = 12 fm. It then starts decreasing with larger b. This is easily understandable. 
There are basically two factors determining the magnetic field strength for a given beam energy. One is the position vector 
R from the moving charges to the field point, and the other one is the charge number of the spectator N^- Their competition 
determines the strength of the magnetic field. For head-on collisions, equivalently there are two counter currents leading to an 
almost zero magnetic field at the center of the reaction. For off-central collisions, as the impact parameter increases, while the 
spectators are farther away from the center they carry more charges. The net result is that the magnetic field becomes stronger 
with increasing impact parameter However, as the impact parameter becomes larger than the sum of the radius of the projectile 
and target, e.g., when b> 12 fm for the Au+Au reaction, almost all charges are with the spectators, the magnetic field is thus only 
determined by the R. Therefore, the reactions with lager impact parameters create weaker magnetic fields at the center of the 
reaction. Based on the IBUUl 1 results, off-central collisions with /? = 8 ^ 10 fm seem to be the most suitable impact parameter 
range to produce the strongest magnetic effect. These reactions create strong magnetic fields and also enough light charged 
particles moving in the magnetic fields to be detected in experiments. Another factor determining the strength of magnetic field 
is the velocity of spectators, i.e., the beam energy of the reaction. Shown in the left panel of Fig. 4 is the beam energy dependence 
of eBy{Q). As one expects, while the maximum strength of the magnetic field increases with beam energy the duration of the 



strong magnetic field decreases since the spectators leave the collision region quickly at higher beam energies. Compared to 
reactions at RHIC, the strength of the magnetic field is about 10 times lower but the reaction lasts about 10 times longer. Since 
observable effects of any force depend on not only its strength but also its duration, magnetic effects in HICs at intermediate 
energies are thus worth an investigation. 

B. Magnetic effects on observables in heavy-ion collisions 

While no chiral magnetic effect is expected in HICs at intermediate energies, it is still interesting to examine magnetic effects 
on hadronic observables. First of all, we would like to mention that the effects of strong magnetic fields on the Equation of State 
(EOS) of cold hadronic and quark matter including the Landau quantization and the nucleon anomalous magnetic moment in 
neutron stars have been studied extensively, see, e.g., refs. [22-26]. It has been shown consistently that the magnetic effects 
become significant only for magnetic fields stronger than about lO"* G. Moreover, at finite temperature some of the magnetic 
effects get mostly washed out [26]. Since the temperature is high and the maximum strength of the magnetic field created is still 
below lO'^ G even at RHIC energies, it is not necessary to consider effects of the magnetic field on the nuclear EOS. Instead, we 
focus directly on magnetic effects due to the Lorentz force acting on moving charges. In the following, we examine separately 
magnetic effects on nucleons and pions. 

1. Lorentz force compared with the Coulomb and nuclear forces 

For nucleons, the magnetic effects are expected to be negligible as the Lorentz force is known to be very small compared to 
the nuclear force. On the other hand, while the electrical and magnetic fields are strongly correlated, the Coulomb force has 
been routinely taken into account but the Lorentz force is normally neglected in modeling HICs. To check the validity of this 
practice and obtain a more quantitative understanding about the relative importance of the Lorentz, Coulomb and nuclear forces, 
we examine in Fig. 5 the ratios of the Lorentz force over the Coulomb and nuclear forces for a test-charge. To be specific, we 
calculate the ratio R^^ of the A:-component of the Lorentz force over that of the Coulomb force for a test-charge in the outside- 
zone. As a reference, we first make an analytical analysis for a simplified case. For a test-charge located at the surface of the 
projectile moving on the trajectory of t(— ^ — R, 0, zo + vqI), where R, zq and vq is the radii, initial z-coordinate and the beam 
velocity, assuming the electromagnetic fields are due to two moving point charges (projectile and target) given by Eqs. (4) and 
(5), the Rf'^ is simply 

/.-^^^f^^f^f. (6) 
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Thus, it is clear that only for fast moving particles likely existing in reactions at high beam energies, the Lorentz force is 
expected to be significant compared to the Coulomb force. We now examine numerically the R'^^ for the test-charge using 
the electromagnetic fields calculated with the IBUUIL In window (a), the time evolution of R'^^ is shown for several impact 
parameters for the 500 AMeV AuH-Au reactions. The evolution can be approximately divided into four periods. Before the 
two nuclei get in touch, R'^^ =0.21 which is exactly the same as the prediction of Eq. 6. In the compression phase, since the 
magnetic fields in the outside region generated by the projectile-like and target-like spectators are in the opposite directions, the 
net magnetic field decreases whereas the electric field there becomes stronger. Consequently, the R'f^ drops until about 15 fm/c. 
In the expansion phase, the situation is reversed. After the collisions are over, the R'^^ keeps approximately a constant value 
smaller than (I'o/c)^ depending on the impact parameter. The beam energy dependence shown in window (b) for the AuH-Au 
reactions with ^=10 fm can be similarly understood. We notice that R^^ = (vq/c)^ at the beginning of the collision is satisfied 
at all beam energies. As the incident energy increases, the Lorentz force becomes closer to the Coulomb force. 

We now turn to the ratio between the .v-components of the nuclear and Lorentz forces, i.e., 7?^*' = F^ /F^, for a test-proton 
at the center of mass with a constant velocity of Vj = vq. Shown in windows (c) and (d) are the impact parameter and beam 
energy dependences of the R^'^ . Because the nuclear force is proportional to the gradient of the single-nucleon potential, i.e., 
F'^ = —VrU, large fluctuations are seen in the R'^'^ . It is seen that the nuclear force is several 10 to 10^ times larger than the 
Lorentz force. The magnetic field is thus not expected to affect the reaction dynamics and nucleon observables. Therefore, it is 
not surprising that nuclear reaction models can describe most experimental data without considering any magnetic effect at all. 

2. Magnetic effects on collective obsen'ables of nucleons and pions 

While the magnetic effects on nucleon observables are expected to be very small, to be quantitative it is still necessary to 
examine how small the effects are. From the expression of the Lorentz force F'^ — qv x B, it is easy to see that the main 
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FIG. 5: (color online) The ratio between the x-components of the magnetic and electric forces (R^^) (a) at various impact parameters for the 
500 AMeV Au+Au reactions, (b) at various incident energies for the Au+Au reactions with b=lO fm; and the ratio between the x-components 
of the magnetic and nuclear forces (Rx'^) (c) at various impact parameters for the 500 AMeV Au+Au reactions, (d) at various incident 
energies for the Au+Au reactions with b=lO fm (the lines are the results smoothed with the Fast Fourier Transformation Filter to guide the 
eye), respectively. 



component of the Lorentz force is in the reaction plane (especially in the jc-direction). The average transverse momentum in 
the reaction plane, i.e., < p^ >, is thus a good candidate. Shown in the top panels of Fig. 6 are the average in-plane transverse 
momentum as a function of rapidity, the so-called in-plane transverse flow [27], for free protons and pions, respectively. Indeed, 
there is essentially no magnetic effect on nucleons. It is seen that both negative and positive pions flow in the same direction as 
nucleons but with much lower transverse momentum in the reaction plane [28]. Interestingly, there is a very weak indication of 
some magnetic effects on the < px{y) > of pions at forward/backward rapidities. This is qualitatively understandable because 
the Lorentz force influences pions motion easily as they are light compared to nucleons. Moreover, it also indicates that the 
magnetic field decreases (increases) very slightly the magnitude of < p^ > for positive (negative) pions at both forward and 
backward rapidities due to the magnetic focusing/diverging effects as we shall discuss in detail in the next subsection. Next, we 
investigate in the lower panels of Fig. 6 the so-called differential elUptic flow as a function of transverse momentum [29, 30], 
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where A^ is the total number of free particles. The piy is /th particle's transverse momentum perpendicular to the reaction plane. 
Again, there is essentially no magnetic effect on the differential elliptical flow of both nucleons and pions. 
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FIG. 6: (Top pane) The average in-plane transverse momentum of free protons and pions as a function of rapidity. (Bottom pane) Elliptic flow 
for free protons and pions as a function of transverse momentum for the 2 AGeV Au+Au reaction at an impact parameter of 5 fm with x=0. 



3. Magnetic effects on the TZ JTl^ and neuteron/proton ratio 



It is well known that the Coulomb force affects significantly the %^ jn'^ ratio in HICs. The so-called Coulomb peak often 
appears near the projectile and/or target rapidities. This phenomenon has been studied extensively both experimentally [31-35] 
and theoretically [37-42] since the 1970's, see. e.g., ref. [43] for a review. However, magnetic effects were not considered 
in any of these studies. While the Lorentz force on pions is normally smaller than the Coulomb force, they have the same 
order of magnitude. Moreover, compared to nucleons pions are light with relatively higher speeds and are thus more easily 
affected by the Lorentz force. Furthermore, there is no nuclear force acting on pions once they are produced at least in most 
model simulations where pions change their momenta only through pion-hadron collisions and the Coulomb field. To our best 
knowledge, theoretical studies on the mean-field (in-medium dispersion relation) for pions are still rather inconclusive [44]. 
Considering all of the above, the magnetic force on pions can be significant. In fact, we expect the Lorentz and Coulomb forces 
to have the opposite effects on the n^ /k^ ratio. Namely, near the projectile/targer rapidity the Coulomb force increase the 
K^ /k^ ratio while the Lorentz force reduces it. Effects of the Lorentz forces on positive and negative pions are illustrated 
in Fig. 7 using the projectile-like spectator as an example. The moving track of the spectator can be regarded as a current. 
Above/below the current, the magnetic field is perpendicular to the reaction plane and points outward/inward. The Lorentz force 
focuses the 7r+ into smaller forward (backward) polar angles while disperses the n^ to larger forward (backward) polar angles. 
So the %^ jn^ ratios at large rapidities are reduced by the Lorentz force. Moreover, due to the magnetic focusing/dispersing 
effect on positive/negative charges, the changes in transverse momentum for particles above and below the current are opposite. 
So the total magnetic effect on the average transverse momentum in the reaction plane is very tiny even for pions. This explains 
why the magnetic effects on the transverse flow < px{y) > and the differential elliptical flow V2{pt) are negligible for both 
nucleons and pions. 

Why is it so important to understand clearly and precisely the electromagnetic effects on the %^ JK^ ratio? One special 
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FIG. 7: A sketch of magnetic effect on the n JTl^ ratio near the projectile rapidity. 



reason is that the %^ jn^ ratio has been predicted as one of the most promising probes of the nuclear symmetry energy at supra- 
saturation densities [45]. While comparisons of transport model predictions [46^8] with exiting data [49] are still inconclusive, 
all models have consistently shown that the n^ JK^ ratio is rather sensitive to the high density behavior of the nuclear symmetry 
energy. The latter is rather poorly known as indicated in Fig. 1. In fact, even the trend of the symmetry energy at supra- 
saturation densities, namely, whether it increases or decreases with increasing density, is still controversial partially because 
of our poor knowledge about the isospin-dependence of strong interaction. To extract reliably accurate information from the 
n^ /n^ ratio about the high-density symmetry energy, it is thus necessary to understand precisely effects from the well-known 
electromagnetic interactions. So, how strong is the magnetic effect on the k^ /n^ ratio in comparison to the symmetry energy 
effect? To answer this question and give a quantitative example, we show in Fig. 8 the 7r^/7r+ ratio as a function of rapidity 
with and without the magnetic field calculated with three different values of the symmetry energy parameter x for the 2 AGeV 
AuH-Au reactions at an impact parameter of b=0 and 5 fm, respectively. In each case considered here, 200,000 IBUUl 1 events 
are used. Comparing the results obtained with and without the magnetic field using any of the x parameter considered, it is seen 
that significant magnetic effects on the K^ /n^ ratio are obvious especially at forward and backward rapidities particularly for 
mid-central collisions. Quantitatively, the K^ /n^ ratio obtained with the magnetic field is significantly lower at forward and 
backward rapidities (polar angles) due to the magnetic focusing/diverging effects on the positive/negative pions as we illustrated 
in Fig. 7. Pions at higher rapidities have larger longitudinal momenta and thus feel stronger Lorentz forces compared to those 
at mid-rapidity. For the head-on collisions, the magnetic effect is small but still appreciable especially in the early phase of the 
reactions when most of the pions are produced. From peripheral to head-on collisions, the ;r^/7r+ ratio changes gradually from 
forward-backward peaked to center-peaked distributions. In peripheral collisions, there are significant Coulomb effects due to 
the spectators. One thus expects the n^ JK^ ratio to peak at forward-backward rapidities. It is seen that the magnetic effect at 
forward-backward rapidities is compatible with the symmetry energy effect from changing the x parameter by one unit. Overall, 
the n^ /k^ ratio decreases as the symmetry energy at supra-saturation densities becomes stiffer when the parameters changes 
from 1 to - 1 . 

It is worth noticing that so far only the integrated 7r^/7r+ ratio, i.e., the ratio of total n^ to 7r+ multiplicities, has been used 
in attempts to constrain the symmetry energy at high densities without considering the magnetic effects. While the integrated 
7r^/7r+ ratio is rather sensitive to the symmetry energy parameters in reactions near the pion production threshold, as the beam 
energy becomes higher than about 1 Gev/nucleon, the sensitivity gradually disappears [46]. It is thus interesting to see that the 
rapidity distribution of the 7r^/7r+ ratio shows a strong sensitivity to the parameter x even in the reactions at a beam energy 
of 2 GeV/nucleon where the baryon density can reach about 3.5p(). Since the strongest sensitivity to the symmetry energy is 
at forward and backward rapidities where the 7r"/7r+ ratio is also strongly affected by the magnetic field, special cares have 
to be taken in both model calculations and the data analysis. Most of the available detectors including the one used by the 
FOPI Collaboration [49] do not provide full coverage at very forward/backward angles. The integrated n^ jn^ ratio is normally 
obtained by extrapolating the angular distributions of pions measured in a limited angular range to all polar angles. By doing so, 
however, the magnetic effects on the angular distribution were neglected. Previous conclusions on the high-density symmetry 
energy based on comparing various transport model calculations with the experimental data without considering the magnetic 
effects thus need to be taken with caution. For comparisons, the neutron/proton ratio njp of free (selected as those with local 
density less than po/8 at freeze-out) and all nucleons are shown as functions of rapidity in the middle and bottom windows 
of Fig. 8, respectively. It is seen that there is essentially no noticeable magnetic effects within error bars on the nj p ratios. 
This is consistent with our expectation and the results on the transverse and elliptical flows discussed earlier The non-uniform 
njp and n^ JW^ ratios as functions of rapidity indicates the lack of complete isospin equilibrium for both the nucleon and pion 



10 



2.1 



2.0 



1.9 



'*^ 1.. 



1.7 




1.4- 



1.3 



0} 
03 



Q- 1.2 



1.1 



1.0 



1.4- 



1.3 



-2- 



1.2 



1.1 



1.0 



x=1 



black b=0 fm 
red b=5 fm 

-■ — without B 
<^ with B 




x=0 



y 

/.^A, 



Au+Au 

E, =2 AGeV 

beam 



\ A 

'\ / 




x=-1 



-1.0 -0.5 0.0 0.5 

y/y, 



FIG. 8: The n^ jlC^ ratio (top windows), the neutron/proton njp ratio of free (middle windows) and all (bottom windows) nucleons as a 
function of rapidity with and without the magnetic field calculated with the three different values of symmetry energy parameter x for the 
reactions of 2 AGeV Au-fAu at an impact parameter b=0 and 5 fm, respectively. 



components. This is the so-called isospin translucency expected in heavy-ion reactions at the beam energies studied here [51]. 

Shown in Table I are the integrated n^ JK^ and neutron/proton ratios calculated without/with the magnetic field. It is seen 
that the integrated ratios are not affected much by the magnetic field. This is what we expected as the Lorentz force affects 
differently only the angular distributions of positively and negatively charged particles, but not their total multiplicities. Also, 
consistent with previous findings [46] the integrated %" JTi^ ratio at beam energies higher than about 1 GeV/nucleon is not so 
sensitive to the variation of the symmetry energy while there is a clear indication that a higher 7r^/7r+ ratio is obtained with a 
softer £'sym(p) at supra-saturation densities. Thus, the differential n^ /n^ ratio as a function of rapidity, as we discussed earlier, 
is a better probe of the symmetry energy at supra-saturation densities after taking care of the magnetic effects. 
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TABLE I: Integrated n^ /;r+ and n/ p ratio calculated without/with the magnetic field using three values of the symmetry energy parameter 
X = 1 , and — 1 . 



Ratio 


b (fm) 


x=\ 


x=0 


A-=-l 


71-/71+ 




5 


imnsi 

1.87/1.86 


1.81/1.78 
1.79/1.79 


1.68/1.67 
1.73/1.73 


n/p (free) 




5 


1.23/1.23 
1.28/1.28 


1.24/1.24 
1.29/1.29 


1.25/1.25 
1.29/1.29 


n/p (all) 




5 


1.23/1.23 
1.31/1.31 


1.24/1.24 
1.31/1.31 


1.25/1.25 
1.32/1.32 



IV. SUMMARY 

In summary, within the transport model IBUUl 1, the time-evolution and space-distribution of internal electromagnetic fields 
in HICs at beam energies between 200 and 2000 MeV/nucleon are studied. While the magnetic field can reach about 7 x lO'^ 
G, it has almost no effect on nucleon observables as the Lorentz force is normally much weaker than the nuclear force. On the 
other hand, the magnetic field has a strong focusing/diverging effect on positive/negative pions at forward/backward rapidities. 
Consequently, the differential 7i^ /n+ ratio as a function of rapidity, but not the integrated one, is significantly altered by the 
magnetic field. At beam energies above about 1 GeV/nucleon, the differential n^ /k+ ratio is more sensitive to the ^symlp) than 
the integrated n^ /n+ ratio. Our findings suggest that magnetic effects should be carefully considered in future studies of using 
the differential 7r^/7r+ ratio as a probe of the E^ymip) at supra-saturation densities. 
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